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Abstract 

We compute the potential between a pair of nucleons in the D4-D8 holographic 
QCD. In the large 't Hooft coupling limit, A S> 1, the hadronic size of the baryon 
is small ~ \j\J XMkk, and their interaction with mesons are well approximated 
by a set of dimension four and five operators. The nucleon-nucleon potential 
emerges from one-boson exchange picture involving massless pseudo-scalars 
and an infinite tower of spin one mesons. We find in particular that p meson 
exchanges are dominated by a dimension five derivative coupling of tensor type, 
whereas for ui mesons and axial mesons, such tensor couplings are completely 
absent. The potential is universally repulsive ~ 1/r 2 at short distance, and 
has the usual long-distance attractive behavior ~ — 1/r 3 along a isosinglet and 
spin triplet channel. Both the large N c form and the finite N c form are given. 
In the former, a shallow classical minimum of depth ~ 0.1MkkN c / A forms at 
around rM^K — 5.5. 
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1 Nucleon-Nucleon Potential from Holography 

The main goal of this paper is to extract the interaction between a pair of nucleons 
in string theoretical framework of holography, and consider possible bound state. 
We will exclusively work with the D4-D8 model [TJ, which involves a large number 
of colors N c , large 't Hooft coupling A, and quenching of fermions. Given many 
approximations, the result should be approached with much caution, yet we have 
seen often that such holographic approaches generating realistic numbers. Ref. [2JE], 
for instance, gave detailed predictions on glueball spectrum of pure QCD, some of 
which were successfully compared to lattice simulation. 

The D4-D8 model has been particularly successful in encoding the spin 1 meson 
sector coupled with pseudo-Goldstone bosons, and equally successful in describing 
baryons and the interaction between these two sectors [U HJ El E] . A natural extension 
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of the model would be a study of the nucleon-nucleon potential. The holographic 
baryon in the D4-D8 model is similar to the Skyrmion of chiral perturbation theory. 
In fact, one can view the holographic baryon as a direct uplift of the Skyrmion in the 
holographic sense, which has a simple interpretation as instanton solitons with certain 
Coulombic electric hair. The size of the soliton is known to be ~ 9.6/ vXMkk [HE], 
where M^k is some natural unit in the D4-D8 model, comparable to the lightest 
vector meson mass. 

Given this, one can approach the problem of nucleon potential in two different 
manners. The first, which seems conceptually most natural, is to find a "suitable" 
family of two-soliton trial configuration, emulating a pair of baryons (nucleons) taken 
apart from each other, and evaluate the resulting energy. After subtracting twice the 
mass of the baryon (nucleon), this would give us a potential. However, this is much 
easier said than done. The main problem with approaches like this is that finding a 
"suitable" configuration is all but impossible for complicated solitons like this. As 
we will see later, the baryon-baryon interaction scales as N c / A whereas their masses 
scale as N C X, so the interaction accounts for a very small part of the two body energy. 
Unless our trial configuration is extremely fine-tuned, energy cost due to any sli ght 
error could easily overrun the interaction energy, resulting in a nonsensical answer 
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Sometimes, however, the interaction energy grows substantially and the physical 
mechanism responsible for the interaction is easy to single out. For our solitonic 
baryon this happens when the two baryons approach each other to a distance com- 
parable to their individual soliton size. The leading contribution comes from the fact 
that each unit soliton comes with Coulombic hair. Each baryon has N c unit of charges 
with the squared electric coupling ~ l/N c X, so one finds a repulsive core interaction 
of type 

v Nc 1 n n 

Vcore ~ TIw (L1) 

with r 2 ~ \jM\ K or less. Here r denotes the mutual separation of the two baryons 
and the 1/r 2 behavior originates from the fact that the soliton lives in approximate 
R i+1 . 

Nevertheless, this short distance behavior gives little insight to some common 
questions like how bound states, such as deuteron and other nuclei, form. In fact, it 
is not clear whether there is a realistic regime where the precise functional behavior 
of the repulsive core such as this can be measured, since, at short distances, the 
asymptotic freedom takes over and nucleons begin to see each other as collection of 
partons. Also the D4-D8 model, or any other holographic model based on gravity 
only, becomes somewhat dubious in the high energy regime well beyond M^k because 
of many non-QCD modes that begin to populate at Mkk and higher. 



#1 In fact, approaches of this kind have been already tried for Skyrmions with mixed results [3 [8]. 
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If one is interested in longer distances where the interaction is potentially attrac- 
tive and where, more to the point, the validity of the present approximation can 



meet real QCD, we must consider a different approaches When the inter-baryon 
distance is larger that the sizes of the baryons, the baryons can be taken to be a 
point-like object. In such situation, the interactions can be all ascribed to exchange 
of light particles, namely mesons. Instead of trying to understand intricate structure 
of multi-solitons, one merely computes Feynman diagrams using (cubic) interaction 
vertices involving baryon currents and light mesons, such as pions tt and rho mesons 
p. Typical vertices that enter this computation are 

rf(x)T(j)(x)Af(x) 

for the (pseudo-) scalar mesons, and 

A/"(x) 7At r^(x)7V(x), N(x)Y u Td^v u (x)N{x) 

for (axial-)vector mesons, with T = 1 or 7 5 . From the D4-D8 holographic QCD, 
coupling constants for these vertices are all precisely derivable, at least in the large A 
and large N c limit. Then the problem of finding nucleon-nucleon potential becomes a 
matter of computing summing up tree-level Feynman diagrams due to various meson 
exchanges 0QISCEI]. 
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Fortunately, the basic framework for the relevant meson-nucleon interaction has 
been worked out in great detail [3], where all nucleon- meson couplings can be de- 
rived very precisely. Some of the leading interaction strengths, such as the leading 
axial coupling to pions and vector meson couplings g p MM an d g^MM, have been 
computed and successfully compared to experimental data. In this note, we will take 
this holographic formulation of nucleon-meson interaction and compute the nucleon- 
nucleon potential from exchange of mesons. For distances r 3> \^X/Mkk, we find the 
leading potential of type 

N c 

^exchange ~ ""T - (' * * ) j (1*^) 

where the ellipsis contains terms of order 1/r through 1/Mj^r 3 possibly with expo- 
nential damping factor due to vector meson masses and also with various spin/isospin 
factors. Subleading contributions start at 1/N C X, but could be relevant in the real 
QCD regime of N c = 3. 

However, this is not to say that the underlying mechanism for the interaction 
energy is different from what one would have obtained from the soliton approach 
if the latter were possible at all. Rather, the tree diagrams involving mesons and 



# 2 The reverse is also true. One should not be tempted to use the second approach for the short- 
distance interaction, since, with much more energy involved, the deformation of individual solitons 
is inevitable. 

# 3 See Ref. [TJ] for a comprehensive review in the conventional approach to QCD. 
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baryons keep track of the classical effect on one baryon by another far away, and 
vice versa, and extract the classical interaction energy automatically. This is possible 
because the baryons are really solitons made out of these mesons, to begin with. 

In section 2 and 3, we give a bare-bone review of D4-D8 holographic QCD and 
the solitonic baryon thereof. In section 4, we derive relevant meson-nucleon couplings 
with emphasis on how they scale with A and N c . In doing so, we will learn that certain 
derivative coupling of vector mesons, sometimes referred to as the tensor coupling, 
can be dominant over the usual minimal type couplings. We extract the values of 
these couplings in the large N c limit. 

Section 5 and 6 summarize the resulting nucleon-nucleon potential, after suitable 
truncation by mass, for large N c and finite N c respectively. It will become clear that 
the leading contributions come from exchange of ir, u, and axial-vector mesons via 
the minimal couplings, and from exchange of p mesons via the tensor coupling, uj 
meson exchange is universally repulsive and represents a remnant of the core repulsion 
we mentioned earlier, while p exchange off-set much of the pion exchange. For very 
short distances ~ /\/XMkk, where the current approach become unreliable due to 
backreaction of the individual solitons, the potential turns universally repulsive as 
~ 1/r 2 as noted above. 

Section 7 gives a simplistic view of deuteron emerging from the large N c form of 
the potential, and we close with a summary in section 8. 

2 A D4-D8 Holographic QCD 

One starts with a stack of D4 branes which is compactified on a thermal circle [13J, 
where one requires anti-periodic boundary condition on all fermions along the circle. 
The purpose of having a spatial "thermal" circle is to give mass to the fermionic 
superpartners and thus break supersymmetry. By putting N c D4 branes on a thermal 
circle, we obtains pure U(N C ) Yang-Mills theory in the remaining noncompact 3 + 1 
dimensions. We are interested in large N c limit, so the U(l) part can be safely 
ignored, and we may pretend that we are studying SU (N c ) theory instead. One then 
extrapolates the AdS/CFT [14j to this non-conformal case, which states that, instead 
of studying strongly coupled large N c Yang-Mills theory, one may look at its dual 
closed string theory. The correct dual geometry is known to be [TS] 

ds 2 = (^Q 3/2 {^dxV + f{U)dr 2 ) + (^ 3/2 Q^l + U'dQ^j , (2.1) 

with R 3 = ng s N c l 3 and f(U) = 1 — Ux K /U 3 . The topology of the spacetime is 
R 3+1 x D x S 4 , with the coordinate r labeling the azimuthal angle of the disk D, 
with t = t + 5t and 8r = AirR 3 ^ 2 /(SU^k)- The circle parameterized by r is the 
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thermal circle. The dilaton is 



while the antisymmetric Ramond-Ramond background field C3 is such that dCs car- 
ries N c unit of flux along S" 4 . 

To add mesons, one introduces Np D8 branes, which share the coordinates 
with the above D4 branes [1] and are transverse to the thermal circle r. If we had 
not traded off the N c D4 branes in favor of the dual gravity theory, this would have 
allowed massless quarks as open strings ending on both the D4 and the D8 branes. 
As the D4's are replaced by the dual geometry, however, the 4-8 open strings have to 
be paired up into 8-8 open strings, the lightest of which belongs to a U (Np) gauge 
field, and these are naturally identified as bi-quark mesons. The U(Np) gauge theory 
on D8 branes has the action 



_ 47^8 J v /Z^ e -* tr^ + ^J C 3 ATre 2 ^, 



(2.3) 



where the contraction is via the induced metric of D8 and \x v = 27r/(2irl s ) p+1 with 
I 2 = a'. The induced metric on the D8 brane is 

r/3/2/ \ D3/2 

hs + x = (dw 2 + n^da?) + jjy^dVl , (2.4) 

after we trade off the holographic (or radial) coordinate U in favor of a conformal 



one w a; 
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w= f R z/2 dU'/Ju' z -Ul K , (2.5) 

JVkk 

which resides in a finite interval of length ~ 0(1/ M KK ) where M KK = 3U^/2R 3/2 . 
Thus, the topology of the D8 worldvolume is R 3+1 x I x S 4 . The nominal Yang-Mills 
coupling g\ M is related to the other parameters as 

9ym = 2ti 9s MkkIs ■ (2.6) 



* 4 This w coordinate is related to another convenient choice of radial coordinate 

U 3 = Uf <K + U KK z 2 , 

as 

d{wM KK ) = ^- = {1+z y ulcKr/3 d(z/U KK ). 
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The low energy parameters of this holographic theory are M^k and A, which together 
with N c set all the physical scales such as the QCD scale and the pion decay constant. 

In the low energy limit, this is reduced to a five-dimensional Yang-Mills theory 
with a Chern-Simons term 

' '' 1 i V = ^tr^ 2 + -^ / u 5 {A), (2.7) 



4 Ja+i e(w) 2 v ^ 24vr 
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where the position-dependent Yang-Mills coupling of this flavor gauge theory is 

1 _ ^Vs* _ XN C U(w) 
e{wf 27r(27r/ s ) 5 108tt 3 KK U KK ' 1 ' 

with V54 the position-dependent volume of S 4 . The Chern-Simons coupling with 
djj*,(A) = trjF 3 arises because J s4 dC^ ~ N c . 

The usual Kaluza-Klein reduction results in an infinite number of vector fields, 
whose action can be derived explicitly as 

jdx 4 C = J dx 4 g tr I X -TfjT^ + m\ n) vV>vM } + •••, (2.9) 

with T^v = d^v^ — d v v^' . When we decomposed U(N F ) into SU (N F ) and U (1), the 
natural gauge generators are normalized as tr T 2 = 1/2, which explains 1/2 in front of 



the kinetic termj^j These fields can be seen as non-zero modes in the decomposition 



of the gauge field, which in the (somewhat illegal but convenient) axial gauge A w = 



is 



A^w) = za M (x)V>(o)(>) + iPn(x) +^2v^\x)ip {n) (w) . (2.10) 



n>l 



The eigenfunctions ip( n ) obey the orthonormality conditions, 



1 



dw ^ (w )(w)*-0(m)(w) = S nm , (2.11) 



2e(w) 2 

for n, m > 1. For later purpose, it is useful to introduce 



/2167T 3 

VV)(™) = \l -jjf-A»)( W ) ( 2 - 12 ) 



# 5 In the published and all prior versions of Ref. [6], the kinetic terms of vector and axial- vector 
mesons were normalized with 1/4 in front of the kinetic term before the trace. With canonical 
normalization for (axial-)vector mesons, the cubic couplings involving a vector or an axial-vector 
meson there should be all multiplied by 
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whose form is insensitive to XN C . Here and in what follows, w = wM^k- Because 
A has a specific parity, the parity of t> n 's are determined by the parity of the eigen- 
functions i)j^(w) along the fifth direction. Since the parity of any one-dimensional 
eigenvalue system alternates, an alternating tower of vector and axial-vector fields 
emerge as the masses mi n ) of the KK modes increase. 

To understand this zero mode part, captured in part by the nonnormalizable 
eigenfunction, xjj^y it is better to give up the axial gauge and consider the Wilson 
line, 

U(x) = e i ^ A(x > w) , (2.13) 

which, as the notation suggests, one identifies with the pion field U(x) = e 2 ™{x)/U_ 
Upon taking a singular gauge transformation back to A w = 0, one finds that it is 
related to a and /3 as 



C/x (.3', 



= {U^ 2 , d.U 1 ' 2 } , 2(3,(x) = [U- 1 / 2 , d.U 1 ' 2 } . (2.14) 



Truncating to this zero mode sector reproduces a Skyrme Lagrangian of pions [16J as 
a dimensional reduction of the five-dimensional Yang-Mills action, 

f dx * (£ tT (u-%U) 2 + - r ± tr[t/- 1 a^ J £/- 1 9 I/ t/] 2 ] , (2.15) 

J \ 4 ile Skyrrae J 

with f 2 = (g 2 YM N c )N c M 2 KK /54^ and l/e 2 Skyrme ~ 61(g 2 YM N c )N c /54n 7 . No other 
quartic terms arise, nor do we find higher order terms in derivative, although we do 
recover the Wess-Zumino-Witten term from the Chern-Simons term To compare 
against actual QCD, we must fix A = gy M N c ~ 17 and M KK ~ 0.94 GeV to fit both 
the pion decay constant f n and the mass of the first vector meson. 



3 Holographic Baryons 

The five-dimensional effective action for the U(Np) gauge field in Eq. (12.71) admits 
solitons which carry a Pontryagin number 

trFAF = k, (3.1) 



with integral k. We denoted by F the non-Abelian part of JF (and similarly later, 
A for the non-Abelian part of A). The smallest unit with k — 1 carries quantum 
numbers of the unit baryon. 

The easiest way to see this identification is to relate it to the Skyrmion [16] of 
chiral perturbation theory, which is the natural object in the large iV c limit [T7| of 
QCD. Recall that both instantons and Skyrmions are labeled by the third homotopy 
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group 7r 3 of a group manifold, which is the integer for any semi-simple Lie group 
manifold G. For the Skyrmion, the winding number shows up in the classification of 
maps 

U{x) : R 3 -> SU(N F = 2) , (3.2) 
while for the instanton it shows up as winding number at infinity, 

A(x,w — > ±00) = ig±(xy dg±(x) , (3.3) 

with 

g-{x) ^g + {x) : R 3 -»• SU{N F ) . (3.4) 

The relationship between the two types of the soliton is immediate [18] once we 
identify 

U{x)=g-{x?g + {x). (3.5) 

Therefore, the instanton soliton in five dimensions is the holographic image of the 
Skyrmions in four dimensions. We will call it the instanton soliton. 

Unlike the usual Yang-Mills theory in flat R A background, the effective action has 
a position-dependent inverse Yang-Mills coupling l/e(w) 2 which is a monotonically 
increasing function of \w\. Since the Pontryagin density contributes to the action as 
multiplied by l/e(w) 2 , this tends to position the soliton near w = and also shrink 
it for the same reason. The F 2 en ergy of a trial configuration with size p can be 
estimated easily in the small p limit, 
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^Pontryagin = 7^ M KK X ( X + ^ M lxP' + ~ , (3-6) 

which clearly shows that the energy from the kinetic term increases with p. This by 
itself would collapse the soliton to a point-like one, making further analysis impossible. 

A second difference comes from the presence of the additional Chern-Simons term 
~ tr^4 A T A J 7 , whereby the Pontryagin density F A F sources some of the gauge 
field A minimally. This electric charge density costs the Coulombic energy 

1 e(0) 2 iV 2 

^coulomb ^ - x I 1 ; c + • • • , (3.7) 

l W7t z p A 

again in the limit of pMxx "C 1 • This Coulombic energy tends to favor larger soliton 
size, which competes against the shrinking force due to Ep ontr yagin- 
The combined energy is minimized at [U EE IE] 

Pbaryon — ~~ 7= > (3-8) 



# 6 The estimate of energy here takes into account the spread of the instanton density D(x l ,w) ~ 
p 4 / (r 2 + w 2 + p 2 ) 4 , but ignores the deviation from the fiat geometry along the four spatial directions. 
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and the classical mass of the stabilized soliton is 



classical f rp _i_ rp \ 

'•"B l k - C/ Pontryagin ~r -'-'Coulomb J 



minimum 



AiV c / */2-3 5 -7r 2 /5 

- 27^ x ( 1+ a ■ (39) 

As was mentioned above, the size pbaryon is significantly smaller than ~ 1/Mkk- We 
have a classical soliton whose size is a lot smaller than the fundamental scale of the 
effective theoryF^l 

For the sake of simplicity, and also because the quarks in this model have no bare 
mass, we will take Np = 2 for the rest of the note. A unit instanton soliton in question 
comes with six collective coordinates. Three correspond to the position in R 3 , and 
three correspond to the gauge angles in SU(Np = 2). If the soliton is small enough 
(pMxK <C 1), there exists approximate symmetries £0(4) = SU(2) + x S77(2)_ at 
w = 0, so the total rotational symmetry of a small solution at origin is SU (Np = 
2) x SU(2) + x SU(2)_. The instanton can be rotated by a conjugate SU{2) action 
as, 

F -> S^FS , (3.10) 

with any 2x2 special unitary matrices S which span S 3 . Then, the quantization 
of the soliton is a matter of finding eigenstates of free and nonrelativistic nonlinear 
sigma-model onto S 3 [201 |2T] . Under the approximate symmetry SU(Np = 2) x 
SU(2) + x S77(2)_, the quantized instantons are in [22] 

(2s + l;2s + l;l) , (3.11) 

while the quantized anti-instantons are in 

(2s + l;l;2s + l) . (3.12) 

Possible values for s are integers and half-integers. However, we are eventually inter- 
ested in iV c = 3, in which case spins and isospins are naturally half- integral. Thus 
we will subsequently consider the case of s = 1/2 states only, which are nucleons. 
Exciting these isospin comes at energy cost. 



# 7 This tendency of the baryonic soliton shrinking to smaller size can be understood as being due 
to the backreaction of vector and axial vector mesons on the conventional Skyrmion 19 . 
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4 Nucleon- Meson Interactions from Holography 



4.1 General Formulation 

The starting point is the five- dimensional effective action of isospin 1/2 baryons. 
With 

o ( -1 \ i ( <n \ 5 ( 1 



we have the following five- dimensional effective action 
d A xdw 



■i&y m D m B - im B (w)BB + 2 ^p^ B 1 mn F mn B 
d 4 xdw—^ trF mn F mn , (4.2) 

where the covariant derivative is defined as D m = d m — i{N c J^ 1 ^ + A m ) with A m 
in the fundamental representation of SU{Np = 2). The position-dependent mass 
m B (w) ~ 4tc 2 / e(w) 2 x (1 + 0(1/ X)) is a very sharp increasing function of \w\, such 
that in the large iV c and large A limit, the baryon wavefunction is effectively localized 
at w = 0. This is the limit where the above effective action is trustworthy. 

The vertex BFB has the coefficient function, about which we only know the central 
value precisely as 

^ Pbaryon _ N c 1 

3e 2 (0) ~VM'M^ ' 1 ' ' 

which shows that this second interaction vertex can be actually dominant over the 
minimal coupling, although it looks subleading in the derivative expansion. As it 
turns out, this term is dominant for cubic vertex processes involving pions or axial 
vector mesons [5]. How to continue this coefficient function to w ^ is unknown. 
However, for all large AiV c estimate of nucleon-meson interaction terms, only this 
central value matters. We chose to use the specific form above for a later convenience 
but it is important to remind ourselves that the precise choice does not matter. 

To obtain interactions between nucleons and mesons, we mode expand B(x fl , w) = 
B + (x^)f + (w) + B_(x^)f-(w) where 7 5 -B± = ±-B± and the profile functions f±(w) 
satisfy 

d w f+(w) +m B (w)f+(w) = m A ff-(w) , 
-d w f-(w) +m B (w)f-(w) = mj^ff+(w), (4.4) 

in the range w G [—w max ,w max \. The 4D Dirac field for the nucleon is then recon- 
structed as 

M = B + + B_ , (4.5) 
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The eigenvalue is the mass of the nucleon mode J\f(x). Approximating mg(w) ~ 
m da 8 sicai^i + ( w m kk ) 2 /3 + ■••), we find 

m M ~mi assical + 0(M KK ), (4.6) 

so for large A and large N c limit, we can take ~ m c ^ assical ~ \N c Mkk/^7ti '. The 
eigenfunctions f±{w) are also normalized as 

dw |.f_H| 2 = l, (4.7) 

Note that there is a 1-1 mapping of eigenmodes with f-(w) = ±f + (—w), where the 
sign choice is tied to the sign choice for m^. Due to the asymmetry under w — > —w, 
f+{w) tends to shift to the positive w side, and the opposite happens for f~(w). In 
this note, we will take the convention where f-(w) = f + (—w). Both can be taken to 
be real. 

Inserting this into the action (14.21) . we find the following structure of the four- 
dimensional nucleon action 

J dx 4 C 4 = J dx 4 (-iN^d^N - im N NN + £ vec tor + £ axia i) , (4.8) 
where we have, schematically, the vector-like cubic couplings 

Sector = -iN^^N - g^M^vf-^M + g^Y^v^AT , (4.9) 

k>\ k>\ 

and the axial cubic couplings to axial mesons, 
£ ax iai = -^TVa^-E^^ (4-10) 

k>l k>l 

For instance, qa is the axial coupling to pions, whose leading cubic coupling to Af 
appears via 

2i „ 2i 



a: 



SU(2) 



{C\ d^} SU(2) = jd^ + ■ ■ ■ = + • • • . (4.11) 



We will ignore quartic couplings involving more than one spin 1 mesons. 

We must recall an important detail which is suppressed in the notation above, 
regarding the differences between the isospin singlet mesons and the triplet mesons. 
These two are packaged into the five- dimensional gauge field A as the trace part and 
the 577(2) part, respectively. For instance, the vector mesons would show up in A as 

-f- 1} = (f !/ 2 W } + P? ,a y. ( 4 - 12 ) 
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where u's and p's are canonically normalized. A crucial point is that the representa- 
tion of A that appears in the baryon effective action is different from this. Instead, 
the vector meson that enters the baryon vertex has the form 

implying the isosinglet has a relative enhancement factor of iV c . The second difference 
can be seen in the fact that only the isotriplets appear in the BFB vertex in five 
dimensions. In any case, each and every cubic coupling above comes in two different 
varieties, ones for isosinglet mesons, such as rj and oo, and those for isotriplet mesons, 
such as Ti and p. 

4.2 Structure of the Cubic Couplings 

All the coupling constants Qy\, g^vdA-> anc ^ 9 A are calculated by suitable wave-function 
overlap integrals involving f± and ^( n )'s [6j|@ Let us consider the general structure. 
Contributions from the minimal coupling, B^^A^B, has the form, 

dw \f ± (w)\ 2 ^ {n) (w). (4.14) 

A^s contribute to dimension four vertices, most notably Af^p^Af, Af^u^Af, and 
their axial vector counterparts. They also contribute to Af , y fJ, 'y 5 d fl 7TAf, although only 
as a subleading contribution. 

Contributions from BFB have the general forms 

dw f P b aryon \ f^y f±{w) ^ n){w) , (4 .1 5 ) 

-w max \ oe\w) j 

for B^ U F^ U B, and 

dwl ^T ) \f±(w)\ 2 d w ^ (n) (w) , (4.16) 

for Bj^F^B. The latter two sets contribute only to the isotriplets. B^s contribute 
to the derivative couplings such as Af^ v d^p v N ' . Cj^'s generate the large N c leading 
contributions to vertices involving isotriplet axial mesons, such as Af'y^^d^irAf and 
the minimal coupling to the axial vectors Af^rfa^N ' . C's also contribute subleading 
pieces to vertices like Af^p^Af. 

# 8 A related but different approach to these couplings was later formulated in Ref. [53] which 
adopted the conventional methods used for Skyrmions. Since both are based on the classical solitons 
quantized over the moduli space, the end results should be equivalent. 
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From these, we have the following cubic couplings for isospin triplet mesons, 



triplet 

9a 


= + 


(k)triplet 

9v 


= ^tk-l + ^2k 


g(k)triplet 


= + A\ k , 


g(k)triplet 




(k)triplet 
9dA 


— Zn 2k > 



(4.17) 

where the mesons in these vertices are in the form such as aj^ U ^ = 2i/ f n x (d fl 7c a T a /2) 
and p^T a /2. As we will see shortly, the second terms for the first three lines are 
subleading in the XN C limit, so that we actually have g™ vlet ~ 4Cq~, gW tn P let ^ Aj fe _ 1 , 
^triplet ^ 2C^.. Of these, gW*" p ' e * — o identically, implying that axial vectors have 
no derivative coupling in our approximation. 

For isospin singlets, B and C contributions are absent so we have 



singlet 
9 A 


= 2A+ 


g(k)singlet 


- A + 

— A 2k- 


g(k)singlet 


- A + 

— A 2k ■ 



(4.18) 

The mesons in these vertices are in the form such as aji^ = 2i/f n x ((N c /2)d fl i] > ) 
and (Ay2)uy 

Remarkably, even before we go into any detail, we have a prediction that all isospin 
singlet vectors and all axial- vectors have no derivative coupling in this approximation. 



4.3 Scaling of the Cubic Couplings 

The key fact that allows us to extract large XN C behavior of cubic couplings is that 
mfs{w) ~ l/e(w) 2 is proportional to XN C . Relative to the mesonic eigenfunctions 
■0(n)> f± becomes more and more concentrated at w — 0. The two wavefunctions 
are slightly off-set from the center by the amount ~ ±1/(M KK \N C ) with the width 
of order ~ 1/(M KK \ / \N C ). This allows us to approximate /| or /+/_ by a delta 
function at origin in wavefunction overlap integrals such as provided that the 
integrand does not vanish near w — 0. 
For example, it is easy to see that 



/2167T 3 - 

^2fc-l "> ^(2fc-l)(0) = \l -jjf- ^(2fc-l)(t& = 0) (4.19) 



v c 
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in the large XN C limit. Here we also used the fact that ^{n){w = wM KK ) for n > 1 
obey 

e 2 (0) - 

dw y i/j^ n) (w)ip (m) (w) = 5 nm (4.20) 

and are independent of AiV c and of M^k- I n particular, numerically we find ^(i) (0) ~ 
0.597. This number is an important ingredient of the low energy nucleon-nucleon 
potential as we will find later. 

^2fc' s > whose integrands vanish at w — 0, take more care. Using reality and the 
eigenmode equation for f±, 

*lk = —f /±H W» +m B (w)U(w)) V( 2fe) H . (4.21) 

J w 

Since ^(2fc) is odd, the leading contribution arises from the derivative piece, and we 
find 

^-^M^ A ^ oh <4 - 22) 

which scales as l/(AiV c ) 3 / 2 . Note that \A 2k -i\ ~ (AA^)" 1 / 2 > \A 2k \ ~ (AiV c )- 3 / 2 . 

Evaluation of -B's is simpler because it involves /+/_ which is an even function, 
so that 



r 2 n 2 

r> ± rbaryon 
D„ = 



3 

With even ^'s, this gives 



/ i/TWV±(%)W- (4-23) 



D ± 2tt 2 p 2 /2167T 3 ? 

- -3^-y^(*-i)( o )> ( 4 - 24 ) 

whereas the odd cases vanish identically 

B& = . (4.25) 

Evaluation of proceeds similarly as A^. Using the equation of motion for f± 
again and recalling that 

2tt 2 o 2 n 2 
3e(wy 6 

we find 

C± = %^ / / ± W ± I%/ T )W W H. (4.27) 
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This gives 



ri± . Pbaryon A/T 2 / 2167T 3 j „. . 

C 2fc-1 -> ± 12 M KK\I -}Jj-V(2k-l) (0) 



C£ - ^^M^J^^)'(O). (4.28) 



Also note that |C 2A; | ~ y/N c /\ > \C 2k -i\ ~ l/\/A 3 iV c . 

The case of n = requires special attention since ^(o) (w) is not normalizable and 
only its derivative, which is normalizable, appears in the physical quantities. The 
conventional choice is such that 9„^(o)(0) = l/n, which is necessary for the familiar 
chiral Lagrangian to emerge from this formulation. With this, we find 

At - ±— - , (4.29) 

and 

B± = . (4.30) 

Finally, with the specific functional form ~ l/e(w) 2 of the &FB coefficient, we have 
an analytical result, 

± = p\ aryon m M M KK N c 1 

° ~ 6tt ~ v^O 7T • 1 ' j 

These enter pion-nucleon couplings, which come with additional factors of 1/ f n for 
each pion. 

Note that some of the above integrals have signs sensitive to the choice of f±. 
Since f± are wavefunctions specific to the chiral and the anti-chiral spinors, these ± 
signs for the values of A 2k and C 2k _ x have the net effective of introducing a 7 5 to the 
vertex as a part of dimensional reduction process, in addition to the existing Dirac 
matrices of the vertices in (I4.2p . This is already manifest in how these coefficients 
contributes to the cubic couplings in Eq. (14 .171) . (14. 181) . 

4.4 Pseudo-scalar Mesons: it and rf 

Starting with 

-^ATfVcvAA, (4.32) 
we restore the isotriplet and the isosinglet mesons and find 

triplet singlet » r 

A/VV<y 7rV Q )Af + 9a (4-33) 
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Since we will be considering Nf = 2, the distinction between i] and 7/ becomes a 
bit ambiguous. Here rj' denotes the trace part of the pseudo-scalar, regardless of the 
number of flavors. In turn, this is equivalent to 



triplet 



\ / singlet *t \ 

9 f . , J AT 7 5 (7rV a )Ar - ( !l \ 2f ' x 2„, A - j fj^'N. (4.34) 

4.5 Vector Mesons: p and u; 

We will denote the isotriplet vectors by and singlets by uj^ k \ upon which 

- E + £ g^ST-Td^- 1 ^ (4.35) 
fe>i fe>i 

separates to 

((k)triplet\ / (k)triplet\ 

— 2 M^ P f a T a N + £ ^ N^d^^N (4.36) 
/ fe>i V ~ / 



and 

^(k)singlet 

~2 



~ E ( —0 ) (4.37) 

fc>i 

since the singlet does not have the derivative coupling in this approximation. 

4.6 Axial Vector Mesons: a and / 

Similarly, the axial vector mesons couplings 

" E gfXl^v^N + E 9 { ^Y u l%v^ k W (4.38) 

k>l k>l 

can be written as 

(k)triplet \ / AT (k)singlet 



.tfyy a (*) _ £ M^f^N (4.39) 

/ fc>i V / 

since no derivative coupling exists for axial vectors in this approximation. 
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5 Large N c Nucleon-Nucleon Potential 

Phenomenologically the nucleon-nucleon (NN) potential is well described by one bo- 
son exchange models. The long-range part of the NN potential is mostly due to 
the pion exchange, while the short-range repulsion is governed by the vector me- 
son exchange. The scalar meson exchange is responsible for the intermediate-range 
of the potential. The interaction Lagrangians for boson-nucleon couplings are, for 
pseudoscalar mesons: 

£p = -g v MMtt(x) n fc<p(x)N(x) , (5.1) 

and for vector mesons: 

C v = -9vNNrt(x)Yv»{x)M{x) + ^Sf( x )^d^{x)M{x) , (5.2) 
where is the nucleon mass. For the D4-D8 holographic model, we saw that the 



derivative coupling is absent for the isospin singlet vectors such as cuj^j The same is 
true of axial vectors, so we have only [11 J 



C-a = -9aAfAf-Af(x)Yl5a^(x)J\f(x) . (5.3) 

Note that we now use the convention for isovector bosons as (p = f ■ 0, v = f ■ v, and 
a = t ■ a. 

It is useful to compare our convention to that of Ericson and Weise [32], which 
is our primary reference on one boson exchange potential. The Dirac matrices we 
used are such that ry M = 7E r i CSOn _weise; w hich brings us to the same convention for 
the nucleon field and its free Lagrangian. In addition, we have reversed the overall 
sign of the couplings from theirs as g^tf = -g P , g p MM = ~9v, and 9 P MM = ~9t, 
which is a matter of a common sign convention on meson fields. We have no scalar 
field, so do not have the counterpart of their g$. 

The leading large N c and large A scaling is such that, for pseudo-scalars {tp = it, rf) 

triplet 

^AW M = 9_A M 

2rrw 2U 




tit- singlet 

9 -^M KK = f^L—MKK^^/^x^r-^, (5.4) 



# 9 Note that, empirically, g/g = 3.7—6.1 for the p-meson (see for example [9]), while for w-mesons 
the ratio is close to zero, for instance g/g = 0.1 ± 0.2 [12] . 
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for vectors (v = p^ k \uj^) 



!k)triplet 

9v 



iy (k)singlet Hw - 



~2m^T KK = 2 -V^^^ fc - 1)(0)X VT' (5 - 5) 

and for axial vectors (a = a^ k \ /^), 

0a«JW = 3 y g V>(2fc) (0) x y — , 

9fwjw = 2 V ~^T~ ^ 2k) (0) x \N c \lT' (5 ' 6) 

Note that and g p j^j^ we have derived from the D4-D8 model are of the same 

sign, which is consistent with experimental results. 



5.1 Holographic Nucleon-Nucleon Potential 

In general, the one-boson exchange nucleon-nucleon potential can be written as 

oo oo oo oo 

v n + + v> + Y V ^ k) + Y v - (k) + Y v f (k) • ( 5 - 7 ) 

k=l k=l k=l k=l 

We now borrow results on one-boson exchange potentials from Ref. [T0l[T2j for various 
mesons, and truncate to the leading contributions in 1/N C and in 1/A. (For more 
complete forms of one boson exchange potential, we refer to Appendix 10 of Ref. |12j.) 
In doing so, we find that not all terms in the above expansion contribute at the leading 
order. The leading contributions arise from the following four classes of couplings 

g^MM^-KK g p (k)MAf M KK [N c 

~ ■ g^WAfM ~ o ~ 9aWjVAf ~ V ~T ' \ b -°> 



whereas g v >NN is further suppressed by l/\N c and g p wj^j^ by l/N c . 
For instance, the one pion exchange potential (OPEP) would be 

ifrrf~J i2^^°( m7rr )^ 1 ' ^ 2 + y^im^S^n ■ r 2 , (5.9) 
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where 5*12 = 3(er ■ f)(cr 2 • r) — &\ ■ a 2 , and 



e~ x ( 3 3 \ e~ x 
y (x) = , y 2 (x) = ! + - + - . (5.10) 



i ■ \ rf ry* 2 / ) • 



However, since we are working in the D4-D8 model where m n = 0, the OPEP sim- 
plifies to 

v „ = _L (£^£) 2 _J_ Sa n . f2 . ,5.11) 

For the isospin singlet vector meson, namely tu^^-mesons, the derivative coupling is 
absent and the leading large XN C contribution is very simple, 

T rholoqraphic 1 / \2 / \ /r io\ 

v,w = (g^WAf) yoKwj- (5.12) 

For p( fc ) which are the isospin triplet vector mesons, the derivative coupling is domi- 
nant over the minimal coupling. This also simplifies the potential quite a bit as 

■trholographic 



1 / 9 p WAfAfM KK \ m p (k) 



4tt V. 2mjv I 3M : 



2 [2j/ (m ( k) r)ffi • o 2 - y 2 (m p (k)r)S 12 (r)] n ■ 7^5.13) 



The contribution to Vju.) due to the minimal coupling g p j^j^ are suppressed by addi- 
tional 1/N C . 

The potential from exchange of isospin singlet axial vectors is suppressed by 
additional l/(AiV c ) 2 while triplet axial-vector mesons contributes [TT] 

■trholographic 

(9aWNN) 2 l-2yo(m a w r )°i ■ °2 + J/2(m aW r)Si 2 (r)] n - r 2 . (5.14) 

Finally note that the meson masses are all of order M KK and m p {k) = m w (k) < m a {k) . 
The vector masses and the axial vector masses alternate as k increases. 

5.2 Behavior at r ~ 1/M KK > l/^/\M KK 

When the distance in question is longer than 1/Mkk, it suffices to consider contri- 
butions from light mesons only, 

p 

■y _ -y holographic , \ r / y holographic . -^holographic . -yholographic\ 
fe=l 
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where the level p is determined by the short distance scale, down to which we are 
interested. For instance, if we are interested in distance down to 1/(3M KK ), p — 10 
would suffice. 

More generally, with large but finite A, the smallest distance where one can still 
trust this one-boson exchange potential is when the distance is comparable to the 
solitonic size of the nucleon at ~ \ j \P\Mkk- Around this scale, the current set-up, 
where one implicitly assumes each of the unit baryon to be intact, breaks down and 
one must begin to consider backreactions systematically. Thus, although the sum can 
formally extend to p = oo, it is in practice more sensible to cut it off at p ~ a/A/10, 
after taking into accounts various order one factors. 

The relevant (large \N C ) pion coupling is 



9-kMN 
2m_v 



M 



KK 



N c 

U3 Vf' 



while for (axial-)vector mesons we parameterize the relevant coupling as 



ik \/-, ~ 9p{k)MM M KK 
V A 2m N 



Nr. 




u xl T' 



(5.16) 



(5.17) 



Coefficients, (k, Xk, are determined by V'(2fc-i)(0) and ^fc) (*-0> we us ^ f nese values 
in the following table 1, together with the masses (in unit of M KK ) of the vector and 
the axial vector mesons. 



k 


m w (fc) = m p (k) 


^(2fc-l)(0) 


Zk 


Ck 


m a ( k ) 




Xk 


1 


0.818 


0.5973 


24.44 


8.925 


1.25 


0.629 


9.40 


2 


1.69 


0.5450 


22.30 


8.143 


2.13 


1.10 


16.4 


3 


2.57 


0.5328 


21.81 


7.961 


3.00 


1.56 


23.3 


4 


3.44 


0.5288 


21.64 


7.901 


3.87 


2.02 


30.1 


5 


4.30 


0.5270 


21.57 


7.874 


4.73 


2.47 


36.9 


6 


5.17 


0.5261 


21.52 


7.860 


5.59 


2.93 


43.8 


7 


6.03 


0.5255 


21.50 


7.852 


6.46 


3.38 


50.5 


8 


6.89 


0.5251 


21.48 


7.846 


7.32 


3.83 


57.3 


9 


7.75 


0.5249 


21.48 


7.843 


8.19 


4.29 


64.1 


10 


8.62 


0.5247 


21.47 


7.840 


9.05 


4.74 


70.9 



Table 1: Numerical results for masses and coupling constants for spin one mesons inter- 
acting with nucleons. 

In figure 1, we display the shape of the large iV c potential with p — 10 for the iso- 
singlet sector with total angular momentum one and total spin one. By superselection 
rules, the spatial angular momentum is a mixture of and 2, and effectively we have 



'12 



7l ' T 2 = —3, <Ti • <T 2 = 1. 



(5.18) 
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Figure 1: A plot of large N c nucleon-nucleon potential, truncated at p = 10, along its 
most attractive channel, namely isospin singlet, a-spin triplet, and even spatial angular 
momentum with S±2 = 2. The horizontal axis is for the distance, rM^K, while the potential 
energy along the vertical axis is in unit of MkkN c /4it\. 

This is the only channel which is attractive at long distance. All other channels are 
repulsive. See section 7 for more discussion. The minimum of the potential is located 
around 5.5 /M^k which is a little larger than one fermi if we adopt Mkk — 0.94GeV. 
Toward r = 0, the potential becomes repulsive very quickly, and this is consistent 
with the expected short distance behavior we will see in next subsection. 

If we continue past toward small r beyond the region of validity (set by the 
integer p), the potential actually turns attractive again very sharply. However this is 
an artifact of cutting off the series at finite number of KK modes, and should not be 
taken seriously. Indeed, this unphysical turnaround can be seen to occur right below 
1/3M KK where we expect the p = 10 formula to become untrustworthy, at least for 
arbitrarily large A. For finite A, however, it turns out that there is a very simple 
remedy of this problem. The unphysical turnaround turns out to be a combined 
effect of the truncation and certain finite A correction that we ignored in section 
4. By choosing an optimal value of p in accordance with A, one can easily restore 
physical sensible short-distance behavior as we explain in next subsection. 

5.3 Coulomb Repulsion at Short Distance and Finite A Cor- 
rections to the Large N c Potential 

When the distance between the pair of nucleon is much smaller than 1/Mkk and 
comparable to \ j \I~\Mkk-, the above expressions must be summed over all mesons. 
When A is sufficiently large, however, it is clear where the leading contribution comes 
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from. The holographic picture of the solitonic baryon involves an instanton soliton 
with a unit Pontryagin number dressed with Abelian electric charge. When the net 
soliton configuration is smaller than the curvature scale of the background holographic 
geometry, 1/Mkk, the instanton part of the soliton will behave like that of ordinary 
instanton on R A with scale invariance. 

This implies that the leading potential energy beyond the rest masses of the two 
cores should come from the five- dimensional electrostatic energy associated with the 
Abelian electric charge. Roughly each soliton has N c unit of electric charges and the 
five dimensional electric coupling scales as 1/ ^XNc, and this gives repulsive potential 

~Tm^' < 5 ' 19 > 

Details of this potential are, however, more complicated. The electric charge density is 
basically the same as the Pontryagin density, so the precise form of the two-instanton 
solution enters the potential. In particular the relative spatial/gauge orientation 
of the two-instanton configuration must enter the potential, predicting a particular 
spin / isospin-dependence. 

Clearly, the precise and quantitative structure of the short-distance potential can- 
not be captured by the our one-boson exchange potential since the underlying formu- 
lation for the latter ignores the core shape of the soliton other than its spin/isospin 
structures, whereas in the short- distance ~ 1/V\M KK the potential energy is of order 
N c and is comparable to the electric part of the soliton energy. In order to compute 
the precise structure of this short distance behavior, one should at least start from the 
full two-instanton solution, available in the literature either via ADHM construction 



or in the form of Jackiw-Nohl-Rebbi (JNR) ansatz [2"1]P 10 I Unfortunately, however, 
this approach is difficult to extend beyond very short distance, since the analog of 
AHDM or JNR is not available in a curved background. 

Independent of this, as a self-consistency check, we wish to understand how the 
sum over the KK tower of mesons end up producing 1/r 2 behavior at short distance. 
The leading short- distance power from individual meson exchange is 1/r 3 . Since KK 
modes sum over such powers (after taking into account the coefficients carefully) 
cannot make a 1/r 2 form, somehow 1/r 3 terms must cancel in the full summation 
over mesons. For instance, pions contributes ~ 71N c /AttX to the coefficient, whereas 
the p and the first a meson contribute ~ — 80iV c /47rA and ~ 57N c /AttX, respectively. 
Continuing this fashion, one can see that the pion contribution is gradually eaten 
away by the alternating contributions from the pairs (p^ k \ o"). However, the sum 
up to p = 10 can be seen to weaken ~ 1/r 3 from the pion exchange only by a 



# 10 While our work was in progress, there appeared two related papers |25[ 126] that share some 
common goal with our work. The latter in particular worked out a precise short-distance form of 
the potential using ADHM construction of two-instanton. 
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Figure 2: This plot shows the large N c potential at short distance where the naive 
large XN C formula combined with the truncation becomes untrustworthy. Without the 
finite A-correction, the truncated potential turns attractive again at a short distance of 
order tMkk ~ 3/p. The figure is for p = 10. 



factor of half, which is not enough for the anticipated cancellation. In numerical 
plot with p = 10, shown in figure 2, this manifests as an unphysical turnaround at 
rM KK ~ 3/10. 

One reason behind this deficiency lies with the leading AiV c estimate we found in 
section 4. While most of estimate there are safe in the large N c limit, the quantities 
are actually correct only up to ~ 1/A corrections. This comes about because 

2 ^ r f * %^ x m B (w) x (\ - v/2 ' 3 ^' 7r2/5 + 0(A- 2 ) j , (5.20) 

implying that our numbers for g a {k)j^-j^ were overestimated and we must adjust 

(\/2 • 3 5 • 7T 2 /5+ \ 
1-- ^ + 0(\~ 2 )\, (5.21) 

if we wish to understand finite A cases, regardless of iV c ^> 1. 

Let us note that the smallest distance for which we can trust the truncation up to 
the p-th pair is around rM KK ~ 3/p. Comparing this distance against the solitonic 
baryon size, below which the effective theory we used does not make much sense to 
begin with, we find that such a truncated potential should be a sensible approximation 
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Figure 3: This plot shows the large iV c potential V( p= io) now with 1/ A-corrected axial vector 
couplings at A = 1100. The unphysical turnaround at ~ 0.3Mkk disappeared completely, 
allowing a smooth transition to the short distance 1/r 2 repulsive core. 

if we choose p ~ y/X/10. For lar ger p, the idea of point-like nucleon fails as far as 
interaction with heavier mesons are concerned, while for smaller p the potential Vu,\ 
fails at distances far larger than the individual baryon size. 

For A ~ 1000, one may thus hope that the choice p = 10 should be an optimal 
one. As we saw above, however, the naive formula for the potential exhibits that the 
potential begins to fails dramatically at rM KK = 3/p, by turning strongly attractive 
again. Although there is no strong inconsistency with this (since the potential failed 
where it is expected to fail), it looks a little suspicious. The point is that although 
A ~ 1000 seems large, the correction (15.211) to the axial vector couplings represents 
roughly more than 2% reduction and cannot be neglected. What one should do is 
to correct g a (k)j^j^ as in (I5.21j) and reconstruct the potential. Indeed, the numerical 
estimate shows an almost complete cancellation of short distance 1/r 3 when we take 
p = 10 for the case of A = 1100. Figure 3 shows the corrected potential in this 
case, where the unphysical turnaround at tMkk ~ 3/10 disappeared. In the case of 
A ~ 1000, at least, the minimal choice for truncation, p = 10, was also effective. Once 
this leaves behind 1/r 2 terms as the leading short- distance behaviors, the appearance 
of 1/r 2 from the latter via KK mode sum, where the vectors and the axial vectors 
contributes with the alternating sign, follows easily. 

For larger A, the reduction of g a {k)j^^- would be smaller, but at the same time we 
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must sum over more mesons in order to make the potential trustworthy down to the 
distance comparable to the soliton core size. One may speculate that the optimal 
choice is again to sum up to p ~ a/A/10. However, the couplings for large k are 
inherently ambiguous since it depends more and more sensitively to, for example, 
the precise functional form of the coefficient of ETB. This is because the necessary 
mesonic wavefunction ^V)'s are more and more widespread, which also makes the 
couplings prone to systematic errors from how the numerical estimate is cut-off far 
away from w — 0. 



6 Nucleon-Nucleon Potential for Realistic A and 

If one wishes to understand real QCD with N c = 3, one must consider a different 
regime. For instance, we neglected V v i on account of the small ratio 



HiUL^J^Jl J_ <t whenAJVc>>1 (eu) 



c 



in the holographic limit. Yet, if we consider N c = 3 and A ~ 17 (determined by 
measured values of f w ), we find the ratio to be about 2 and is hardly ignorable. The 
estimates here themselves are no longer reliable since we used large XN C limit, but 
this comparison clearly shows us that we cannot expect any small parameter. This 
is in fact a generic problem in going to the realistic limit. 

In computing Feynman diagrams and extracting nonrelativistic potential, another 
small parameter is p/mjy where p is the spatial momentum of the meson being 
exchanged. However, when translated to real space, this ratio can show up either 
as m/m^f or as which is problematic when the meson mass m exceeds the 

nucleon mass. Thus, contribution from exchange of heavy mesons cannot be included 
reliably, forcing us to cut down to pions, r]', p, and uj. Thanks to the universal 
suppression ~ e~ mr for heavy meson processes, this is a good approximation as long 
as we are interested in distances strictly larger that 1/M KK . 

The relevant Nucleon-Nucleon potential is then 

V = V w + V rf + V pW +V uW , (6.2) 

where individual term must be computed as a series expansion of mjm^. Actually, 
the exchange of vector mesons generates a correction to the kinetic term as well, 
the two-body Hamiltonian for a pair of nucleons contains the relative part of the 
Hamiltonian, 

H = - ( — + A ] V 2 + V (6.3) 
Wat / 
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with 



A 



3m ( i) 
16tt 



sLwm + alaw^n • f 2 ) f ^) yo{m m r) , (6.4) 



where 77i(i) = m u (i) = m (i). In the attractive channel, ri • r 2 = —3, as we will see 
later, {g^wj^j^/g^mj^j^) 2 — 14, so the effective reduced mass of this two body system 
becomes smaller as the distance becomes small. 

Let us turn to the potential. V n was already given in Eq. ( 15. lip , while others can 
be inferred from Ref. [12]. The contribution from the trace part is essentially the 
same as the massive pion case (15.91) except the SU(2) generators f\ ■ f 2 are absent 



V, 



i ( g V ' 



ATM 



4tt \ 2mjv 



M, 



KK 



2 rn^i rrirf 



fy (nv0^i ' ^ 2 + y2( m v' r ) S ^] ■ (6-5) 



The mass of r/ is generated by the U(l) axial anomaly, and was computed by Sakai 
and Sugimoto, 

\M KK [Np 



V p (i) is considerably more involved than V (1 
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with 



S = -{&1 + 02) 



Q12 = ~ (Vi ■ L)(a 2 ■ L) + (a 2 ■ L){S X ■ L)) , (6i 



2 

and the spatial angular momentum L. Finally, V^i) is essentially of the same form 
as V (1), except that g u (i) j^jy = and t\ ■ f 2 is absent, 

m (i) 2 if, 1 m (i) 



2?/ (m (1) r)ai • <? 2 - V2{m {1 )r)S 12 (r) 

12 \ m// / 

3 /m ( i)\ 2 yi(m (1) r) -. 1 /m (1) \ 4 y 2 (m (1) r) 



2 J m(!)r 16 \ m A/" J m^r z J 

These are the complete expressions up to the quartic order in terms of spatial mo- 
menta of individual nucleons. 

For N c = 3 and A ~ 17, we found the following numbers that determine the 
couplings here, 

7 04 

4(7+ ~ 0.697, 2A+ ~ 0.136, ~ 5.93, 25+ ~ — — , 2C+ ~ -1.22 (6.10) 
and 

f n ~ 0.0975M^ K , m ( i) = m pC i) = m w <i) ~ 0.818M^, m n , ~ 0.85M^. (6.11) 

The mass m^- has an inherent ambiguity since it would be additively renormalized 
by massive excitations around the soliton. Our definition of the nucleon mass kept 
only one such massive mode, namely the position along w-direction, and according 
to this prescription, we find 

m M ~ l.93M KK . (6.12) 

Unfortunately, the scale of Mkk that fits the physical nucleon mass is about ~ 
500 MeV, as opposed to the one needed to fit the physical p meson mass at ~ 940 MeV. 
This discrepancy between the mesonic and the baryonic scales was previously ob- 
served both in the D4-D8 modelin a slightly different comparison [5] and also in the 
so-called bottom- up approach [28J, and appears unavoidable in the gravity approx- 
imation to the bulk side. We will proceed with these numbers, nevertheless. The 
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Figure 4: A plot of the nucleon-nucleon potential, with N C X = 50 and N c = 3. We again 
drew the potential for the isospin singlet, cr-spin triplet, and S±2 = 2 eigensector, although 
for this finite XN C case we should expect different S12 eigensectors to mix in. This plot is 
only for the purpose of illustrating the general trend. The horizontal axis is rM^K, while 
the vertical potential energy is in unit of Mkk/^- 



couplings that enter the above potential ar a^ 11 



9 -^LM KK = AC °+ 2A °M KK ^4.27, 
2m N KK 2f n 

9-n'NN 2A • N c 

M KK = — — M KK ~ 4.18, 

2m A f 2}n 

+ 2C+ 

9 p wmm = g ~ 2 ' 36 ' 

At -N c 
9uWMM = 2 ~ 8 ' 90 ' 

~ m ™M KK = 2Bt " MkK * 7.04. (6.13) 



2m,M 

Detailed study of this case will be reported elsewhere 



# n It has been observed previously that next subleading correction of some of the operators may 
involve the simple shift N c — > N c + 2 in the leading expressions. This, for example, allows a very good 
match of g^wx with experiment. The origin of this shift, originally suggested by the constituent 
quark models, is not clear from this approach. Here, we chose not to implement this shift but readers 
should be aware that terms from BFB may be affected, leading to quantitatively different numbers. 
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7 Holographic Deuteron: Large N c Results 

In this final section, we explore some basic aspects of deuteron physics with the NN 
potential 

10 

yholographic , \ ^ / yholographic i yholographic . yholographic\ (71) 
7T / j \ p( fc ) Lj(k) a (k) J \ ' J 

k=l 

in the large A and N c limit. To distinguish this from physical deuteron, we refer to 
them as holographic deuterons. 

For a bound state, we need to focus on the long distance attractive channel. The 
large N c potential has a simple spin and flavor structure as 

yholographic = ^ + (ya^ . ^ + ySg^ ~ . ^ ^ 

The massless pion exchange, the dominant contribution in long distance, contributes 
only to Vj? and positively, so an attractive channel requires Sufi • fi < 0. Using the 
fact that S12 acting on cr-spin singlet vanishes identically, and that the nucleons are 
fermions, this forces the isospin singlet (j\ • f 2 = —3) and the a— triplet (&\ • <t 2 = 1) 
channel with even spatial angular momentum. The lowest total angular momentum 
possible is then J — 1, and the positive Su eigensector has the following spatial 
angular momentum mix as 

|L = 1) + V / 2|L = 0)_ 



V3 

In this eigensector, Su = 2. The figure 1 is the plot of the potential in this sector, 
whose classical minimum occurs at 

5 53 iV 

V(r mm ) ~ -0.09AAM KK ^. (7.4) 



■ unit, — ti j- i V iiiiiL / — " ■ " w ^ ^-^n a » 

M KK A 

Note that the binding potential is very shallow. Recall that in the large AiV c limit, 
Mkk is most conveniently determined by the vector meson scale to be around 0.94 GeV. 
Among various scales that enter the baryon energy, we have the hierarchies, 

, ■ , iV 
vim ~ mf™' ~ XN C » E Coulomb ~N c » \V(r min )\~^. (7.5) 

The middle measures the energy related to the classical deformation of the individual 
baryon away from the self-dual soliton, while the last measures the binding energy of 
the nuclei. 

T he fa ct that the nuclei biding energy is small, which is also borne out in real 



world F 12 ! is interesting from the standpoint of the holographic QCD as well. The 



# 12 The physical deuteron has a binding energy of 2.2 MeV [29], which is about 0.12% of its rest 
mass. 
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stringy picture of the baryon says that the individual baryon can be viewed as a D4 
brane wrapped on the compact S 4 of the dual geometry [23 [T]. What we computed 
here is essentially the potential between two such objects separated along the non- 
compact R 3 . The binding energy is positive but suppressed relative to the individual 
rest mass by 1/A 2 , indicating very weak interactions. In terms of the warped string 
scale, a' warped , this power is equals (a' warped M^ K ) 2 . Although the significance of this 
particular power is unclear to us, it does show that the two wrapped D4 branes 
are almost non-interacting at long distances. This seems to suggest that the object 
underlying baryons may remain close to its original BPS nature, despite the super- 
symmetry breaking background of scale M^k and high mass ~ XN c Mkk, which is 
well beyond the cut-off scale Mkk, and may eventually explain why such a high mass 
object is well-described by this D4-D8 holographic QCD. 



8 Concluding Remarks 

In this work, we computed the nucleon-nucleon potential in the D4-D8 holographic 
QCD, which is generated by exchange of five-dimensional flavor gauge field. In four 
dimensional picture, this amounts to exchange of massless pseudo-scalars and an 
infinite tower of spin one mesons. In the large XN C limit, it is sensible to sum up to 
first ~ ^/A/10 vector and axial vector meson pairs, although one may choose to cut it 
shorter according to the shortest distance scale interested. This prescription also gives 
whereto glue the repulsive short distance regime to the more complicated intermediate 
and long distance regime. Some rudimentary aspects of deuteronic bound state is 
explored for large N c case. Consideration of deuteron for realistic QCD regime will 
be explored elsewhere. 

We hope this work will provide a more practical starting point for exploration 
of how holographic QCD fares against experimental data, part of which comes from 
nucleon-nucleon scattering amplitudes. Admittedly, this would involved huge extrap- 
olation to N c = 3 and A ~ 17, where the holographic approach is hardly justifiable 
by the first principle. But, in the absence of any other honest derivation of nucleon- 
nucleon potential, our result should be at least tested against data. In this work, we 
did not attempt to analyze realistic QCD regime and concentrated mostly on large 
N c limit. We wish to come back later to the N c = 3 potential of section 6, and explore 
its consequences. 

Another important application of this work would be in understanding dense 
matter system, such as neutron stars, where the correct equation of state is of some 
importance. In dealing with such a dense system from the holographic QCD, baryon 
density itself were often treated as external input in the form of delta-function density 
in five-dimensions. We hope our nucleon-nucleon potential would allow a more refined 
approach. 
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